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ABSTRACT 

ON" 

On 

Gamma-Ray Burst (GRB) afterglows are commonly interpreted as synchrotron 
emission from a relativistic blast wave produced by a point explosion in an ambient 
medium, plausibly the interstellar medium of galaxies. We calculate the amplitude 
of flux fluctuations in the lightcurve of afterglows due to inhomogeneities in the 
surrounding medium. Such inhomogeneities are an inevitable consequence of 
interstellar turbulence, but could also be generated by variability and anisotropy in a 
precursor wind from the GRB progenitor. Detection of their properties could provide 

£»»«. ■ important clues about the environments of GRB sources. We apply our calculations 

to GRB990510, where an rms scatter of ~ 2% was observed for the optical flux 
fluctuations on the 0.1-2 hour timescale during the first day of the afterglow, consistent 
with it being entirely due to photometric noise (Stanek et al. 1999). The resulting 

q^ [ upper limits on the density fluctuations on scales of ~ 20-200 AU around the source 

of GRB990510, are lower than the inferred fluctuation amplitude on similar scales in 
Q-f the Galactic interstellar medium. Hourly monitoring of future optical afterglows might 

therefore reveal fractional flux fluctuations at the level of a few percent. 
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Subject headings: Gamma Ray Bursts 



1. INTRODUCTION 

Almost all well-localized Gamma-Ray Burst (GRB) sources have shown afterglow emission 
in X-rays for several hours (see Piran 1999 for a review). Often the emission persists on longer 
timescales at lower photon energies, peaking in the optical on a timescale of days and in the radio 
on a timescale of weeks or longer. This long- lasting afterglow emission is most naturally explained 
as synchrotron emission from a relativistic blast wave, produced by the GRB explosion in an 
external medium (see e.g., Paczyhski & Rhoads 1993, Katz 1994, Meszaros & Rees 1993, 1997, 
Waxman 1997a, b). The emission frequency declines with time due to the deceleration of the shock 
wave (Blandford & McKee 1976) and the corresponding reduction in the characteristic electron 
energy and magnetic field amplitude behind the shock. The external medium could be either the 
interstellar medium (ISM) of the host galaxy (Waxman 1997a) or a precurser wind from the GRB 
progenitor (Chevalier & Li 1999). 



Previous theoretical models of the afterglow emission have assumed for simplicity that the 
density profile of the external medium is smooth, i.e. uniform in the case of the ISM or power-law 
with radius in the case of a progenitor wind. However, this simplifying assumption is not expected 
to hold in realistic situations. The ISM is known to exhibit inhomogeneities due to turbulence, 
and stellar winds may vary in time and in solid angle. Since the emitted afterglow flux depends on 
the instanteneous number of shocked electrons, any density inhomogenties are expected to induce 
temporal fluctuations in the afterglow lightcurve. Observations of these fluctuations could provide 
additional constraints on the nature of the surrounding medium and the GRB progenitor. 

In this paper we derive the relation between the spatial power-spectrum of density fluctuations 
in the ambient medium surrounding the GRB source and the Fourier transform of temporal 
fluctuations in the afterglow flux. We focus on small-amplitude (linear) inhomogeneities, as those 
provide the minimal source of afterglow flux fluctuations. §2 presents the formalism used in our 
derivation, and §3 describes our numerical results. We summarize our main conclusions in §4. 
For simplicity, we consider the case where the unperturbed ambient medium is uniform, as for a 
background ISM. The particle density inferred from detailed modeling of some GRB afterglows is 
in the range ~ 0.1-1 cm -3 (e.g. Wijers & Galama 1998; Waxman 1997a, b), as expected for the 
ISM of their host galaxies. 



2. FORMALISM 

The relativistic fireball produced by a GRB explosion starts to decelerate at the radius where 
the energy of the ambient gas swept by the fireball is comparable to its energy output (e.g. 
Waxman 1997a,b), 

o ir>16 ( -^52*10 \ /-.s 

r c = 2 x 10 I I cm, (1) 

where £52 is the hydrodynamic energy release in units of 10 52 ergs, tio is the duration of the GRB 
in units of 10 s, and n\ is the proton density of the ambient medium in units of 1 cm~ 3 . At larger 
radii, the shock follows the adiabatic self-similar evolution (Blandford & McKee 1976) whereby its 
Lorentz factor declines with radius as 

(r\ " 3/2 

r-r.(-) , (2) 



and 



-E52 



1/8 



r c = 270 -srH • ( 3 ) 

The global spectral characteristics of GRB afterglows are naturally explained in terms of 
synchrotron emission by shock-accelerated electrons from this decelerating relativistic shock (see 
e.g., Wijers, Rees &; Meszaros 1997; Waxman 1997a, b; Sari, Piran & Narayan 1998). 



Figure 1 illustrates the geometry of the emission from an infinitesimal volume element in 
spherical coordinates dV = r 2 drdfid(j) at a distance D from the observer. Here fi = cos 9 and 9 = 
along the line-of-sight. We define the emission coefficient j , to be the power emitted per unit 
frequency per unit volume per steradian in the rest frame of the outflowing material. We use prime 
to denote quantities in the local rest frame of the emitting material, while unprimed quantities are 
measured in the rest frame of the ISM. Note that jy/v 2 is Lorentz invariant (Rybicki & Lightman 
1979). For a spherically-symmetric expansion of material which emits isotropically in its local rest 
frame, we have v = v^/{l — /3/i) and j , = P (y ,r,t)/4iT, where 7 and f3c are the Lorentz factor 
and the velocity of the emitting matter, respectively. A photon emitted at time t and place r in 
the ISM frame will reach the detector at a time T given by 

- T t-'A (4) 



1 + z c 

where z is the cosmological redshift of the GRB and T is chosen such that a photon emitted at 
the origin at t = will arrive at the detector at T = 0. Thus we have (Granot, Piran & Sari 1999) 

P , rj* l + z f 2n ^f\ r 2, P'(^(l-(3fi),r,T z +rfi/c) 

F ^ T) = ^dlJo ^J-^Jo "^ 7 2 (1~/W ' (5) 

where di is the luminosity distance to the GRB, and 7, /3, // are evaluated at the time t implied 
by equation (||). 

Most of the shocked material is concentrated in a thin shell behind the shock front. The 
characteristic thickness of the shell is A ~ .R/IO7 2 in the ISM frame, where 7 = T/y/2 is the 
Lorentz factor of the material just behind the shock. In the following, we will assume that the 
observed radiation originates from the thin shell of thickness A = rjR/j 2 behind the shock, inside 
of which the Lorentz factor, the particle density, and the energy density of shocked ISM obtain 
the following values, 

n = Ajn, e = 47 nm p c , (6) 

where n is the number density of the unshocked ambient ISM in its local rest frame, and m p is 
the proton mass. The actual value of rj depends not only on the hydrodynamics but also on the 
behaviour of the magnetic field and the shock-accelerated electrons which determine the local 
emissivity. 

The volume integration expressed in equation (||) should be taken over the region occupied 
by the emitting shell at a given observed time, as illustrated in Figure 2. Because of relativistic 
beaming, the observed radiation originates from a small angle along the line of sight, 9 < I/7. 
Hence, we can set the upper limit in the integration over 9 to be 2/75, where 75 is the Lorentz 
factor of the shell at point b on the plot. For a given observed time T, the outer boundary abc of 
the integration region is defined by the relation (Granot et al. 1999) 

cT z 
Routcr = 1 - 11 + 1/(16 7 2 ) ' (?) 



where 7 = 7 (.Router)- Photons originating at this boundary are emitted from the front of the shell 
and arrive at the detector at the same time T (although they are emitted at different times). 
Similarly, the inner boundary of the integration region is described by 

cT z 

-dinner = ~ —. j 1 \ / o ' V°J 

1 - H + (?7 + I g)/7 2 

with 7 = 7(-Rinner)> and where the associated photons are emitted from the back of the shell. Note 
that although the emitting shell is thin, the region contributing to the observed flux at a given 
time is rather extended; this has important consequences with respect to the spatial scale over 
which density fluctuations could affect the observed flux. Furthermore, consider a point g = (r, 9) 
inside the integration region. A photon emitted from this point at time t = T z + rfi/c will arrive 
at the detector at time T. Since the emitting shell is very thin, the point g is very close to the 
shock front at time t. Thus the radius of the shock at time t is approximately equal to r. This 
fact will be used in our calculation. 

Next, we derive the local emissivity due to synchrotron radiation. We assume that the 
energy densities of the shock-accelerated electrons and the magnetic field are fixed fractions of the 
internal energy density behind the shock front, e e = e e e , e B = e_ge , and that the shock produces 
a power law distribution of accelerated electrons with a number density per Lorentz factor of 
./V(7 e ) = Kj~ p for 7 e > 7 min , where 



P-^ \ e e e ^ _ (m -i^'„,p-i 

p — 1/ n'm e c 2 



7min = ( ~ 7 ) Z7ZTZ 2 ' K = (P ~ l ) n Tmin • ( 9 ) 



Then the emissivity can be approximated as 

P' = tf l7 /3V/V) 1/3 , «/<*W», (10) 

p' = ff27 (3p+l)/2 /3 2 n (p+5)/4 (l/ ' ) -(p-l)/2 ) J > ^. n) (n) 

where ^ min = i / svn (7min) is the synchrotron frequency of an electron with the minimal Lorentz 
factor, and 

ff - go A y (g ~ ^^ n € B 3 m e c3 no) 

Fl " ^^ X (^-1)0,-2)2/3^^3^7^73' (12) 

(n O^P- 1 , n„ m ( 5 P- 3 )/ 4 „3 

fr,-qqn xl q(p-i)/2 x IP ~ 2 J (p+1)/4 p-1 m P c 

ii 2 - cjy.y x i.y x ( ^ _ ^ _ 1)p _ 2 ^re B e e ( 3p_ 5)/2 (3 _ p)/2 U^ 



'mm' 



are constants. Also note that the peak frequency at which the observed flux peaks is fp ea k ~ ^n 

Subtituting equations (|l0|) and (|il|) into equation (J5|) respectively, and making use of the fact 
that for a highly relativistic system, 9 < 2/7^ <C 1, and so 



1 /i~> 

2V 2 ' 2 



d^^-9d9, i_^«_- + _ (14) 



we get the afterglow flux at a frequency v, 
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F ( U ' T ) = T32 / # / «W / dr — — , v < u peak , (15) 

4ird z L Jo Jo Jn innei 1 2/3fj^ + E.y / 



All emission frequencies under consideration are assumed to be below the cooling frequency, v c 
(i.e. not affected by the rapid cooling of the high-energy tail of the electron distribution), and also 
well above the synchrotron self-absorption frequency, f a (see Sari et al. 1998, for more details). 

At point b in Figure 2, the radius of the shock front is Rb, the Lorentz factor of the shell 

behind the shock is 7&, and the time T when photons emitted at point b reach the observer is 

related to these quantities by 

T = J^_ (17) 

(1 + z) 16 76 V v ; 

Note that this time is different from the observation time, T b s , which is defined to be the arrival 

time of most photons emitted from the shell of radius Rb- This is because most photons are 

emitted from a cone of opening angle ~ I/7& around the line of sight and they suffer a longer time 

delay ~ Rb/2j 2 c than the photons emitted on the line of sight (Waxman 1997c). The observation 

time is 

J^ „ J|_ = 8--1— . (18) 

(l + z) 2 7 2 c (1 + z) K ! 

We can use Rb to normalize equations ( |15| ) and ([16]). Based on the scaling, 7 oc r~ 3 ' 2 , we get 
Using equations (|l7|) and (|l9l), we can rewrite equations (0) and (|8|) as follows, 

Pouter + 87^ 2 ^outcr "1 = 0, (20) 

(16?? + l)X ; 4 niier + 8 7 2 fl 2 X inner -1 = 0, (21) 

where x = r /Rb, X outcr = -Router / 'Rb, an d dinner = Rmnei/Rb- In these notations, equations (|l5|) 
and ([l^) obtains the form 

[l + z)H 1 u 1 / 3 Rl ^ AA f 2hb M [ x °^ <^ n^ 3 ex 3 dx 

^? + 2 ) 



F(y,T) = ± ' \ b dcj> dO u <u peak , (22) 






„, ^ N {l + z)H 2 v-^- i y 2 Rhl- 1 [ 2 * 71 /-a/T* ,„ fXoutcr(e, Jb ) n^+^^Odx 

F(y,T) = -—J °— / #/ d0 / x fo+3V2 ' ^>^ P eak, 

4vrd| 7o Jo Jx inncT (e,j b ) x ( 3p - 7 )/2 ^ + ^ V P+ V 

(23) 



where X- mner (9, r -fb) and X outcv (6,jb) can be obtained by solving equations fl20| ) and (|2l|), and Rb 
and 76 are functions of the observed time T. At a later time T = T + r, the shock front moves 
from b to b, while the Lorentz factor and radius of the shell change to % and Rb, respectively, with 



R b \f) ' 7fe VT, 



(24) 



Equations (|22| ) and (^) can be written in the generalized form 

F(T) = f G(r)n y (r)dr, (25) 

where the density n{r) may fluctuate, and 



4/3 for v < i/ peak 

(p + 5)/4 for i/ > i/ peak . 



(26) 



Thus we get, 

(F{T)F(T + r)) = ( / G(r)n y (f)dr I G(r)n y (r)dr) 

dr f d?G{r)G(f f )(n y (r)n y {f f )). (27) 



The angular brackets in the above equation reflect an average over an ensemble of afterglows 
with identical source properties, exploding at different places in the ISM. The ergodic assumption 
implies that an average over many such systems would be equivalent to an average over time for 
the explosion if the ISM is in a stationary statistical state (Reif 1965). For the ensemble average 
of the right-hand-side, we have 

(n y (r)n y (f)) = (n) 2y ([l + 6(r)] y [l + 5{F)] y ) 

- (nf y (l + V -^p±5\r) + y -^ Y ^6\f) + y 2 5{r)5{r)) 

= (n) 2y \l + y (y - 1)£ + y 2 £(f - fj\ , (28) 



where 



'« = ^r 1 ' < 29 > 



^(r — r) = (5(f) 6(f)) is the ensemble-averaged autocorrelation function of the density fluctuations, 
and £q = (5 2 (r)) = (5 2 (r)). As mentioned before, we consider only small (linear) density 
fluctuations. In equation (|28|) , we have implicitely assumed 5(r) <C 1. 



Similarly we have 

(F(T))(F(T + t))= I dr I drG(f)G(r)(n y (r))(n y (r)), (30) 



( n y(r))( n y(?)) « (nfy [i + y ( y - i)&] (31) 

The statistical properties of the ambient gas inhomogeneities in the vicinity of GRB sources 
are highly uncertain, and so we adopt the minimal number of free parameters to describe the 
autocorrelation function, namely we write 

£(r)=£ exp(--), (32) 

where tq is the scale length of the density autocorrelation function. For simplicity, we ignore 
deviations of the expanding shell from spherical symmetry. Our calculation focuses on scales 
much smaller than the size of the emission region (see §3) and so the cummulative effect of many 
small-scale patches of density perturbations averages out during the expansion history of the 
shell. At different points inside the shell, the particle density and energy density [described by 
equation (g)] might be temporarily above or below their average values, but the total energy is 
conserved. 

We can now define the autocorrelation function of the temporal fluctuations in the afterglow 
flux as 

CM - (A(T)A(T + r)> - { f^ T \ T l ~ 1. (33) 

where we consider fluctuations on timescales much shorter than the evolution time of the afterglow 
lightcurve, t <^T, and where 

_ fW-OT) (34) 



(F(T)) 



We now define 



^ 3 fe + TJ for ^<*Vak 

g(x,9, 7b ) = { ^ Jb / N - (p+ 3)/2 (35) 

0x -(%-7)/2 ^ + *j for u>Upcak . 

One of the integrals of interest is 

r2iT r^ht rX ou t C r (0,y b ) 

I\ = / d(j) I d8 I g{x, 0, ~fb)dx 

Jo Jo Jx inncl (0,"fb) 

p2n „ /•2/Ti, „ /•-X'outer (§,%) 

x dxj> I dd g(x,e,j b )dx, (36) 

Jo Jo Jx innci (e,7 b ) 

where X mner (9 ,%) and X ou t er (8,%) are found by solving equations ( p0[ ) and (^). A second 
relevant integral is 

d<p d9 g(x,6,^b)dx 

o Jo Jx innci (e,-y b ) 

/• 27r _ /•2/7b - fX outer (0,7b) 

X d(j> del _ y 2 £(\r-r\)g(x,9,%)dx, (37) 

Jo Jo Jx lnneI (0,%) 



where 

\r — r\ = [(xRb sin 9 — xRb sin 9 cos (f>) + (xRb sin 9 sin 0) 

+(xi? fe cos6»-x^ fe cos^) 2 ] 1/2 . (38) 

Using the above integrals we may write 

C(r) = (A(T)A(T + r)) = / 2 // 1 . (39) 

3. NUMERICAL RESULTS 



We evaluated numerically the integrals in equations (|36|) and ( |37D using a Gaussian quadrature 
method. For each of the two different frequency regions {y < f pea k and v > f pea k)> we considered 
six different cases with values of jb of 100, 22 and 3, which correspond to peak afterglow emission 
in the X-ray, optical and radio wavelength regimes, and for each 7^, we considered two different 
cases of r^/Rb = 10 -3 and tq/ Rb = 10~ 2 . The values of 7&, Rb and f pe ak) the frequency at which 
the observed flux peaks at a given observation time T b s , are given by (Granot et al. 1999), 

-^52\ /^obs,days\ 



7b = 7 . 96 (^£j (fiqpj 1 ) ( 4 °) 



R b = 3.29 x 10 



17 



E^oT, 



52 J obs,days 



1/4 

cm, (41) 



0peak(p) f(p) „l/2,2 pl/2^-3/2 
>peak(2.5)/(2.5)' 

where <^ pea k is a slowly decreasing function of p, /(p) = [{p — 2)/(p — l)] 2 , and Todays is the 
observation time in days and is related to T by equation (|i~8|). 



_7in„ in!5, /T— : ypcakW j W 1/2 2 ^1/2^-3/2 „ ,,, 9 n 

pe ak - 7.29 x 10 VI + z- — — e B e e h b2 -/ obS]days Hz. (42) 



In our numerical calculations we assume p = 2.5 (Granot et al. 1999, Sari, Piran & Halpern 
1999), €b = 0.1, e e = 0.1 and 77 = 0.1. Table 1 shows the associated T, T b s , Rb arid ^ pea k for the 
three choices of 7&, with .E52, ni and 2 as free parameters. 

Figures 3 and 4 show the numerical results for the six cases mentioned above, with Figure 3 

corresponding to v < z/ poa k and Figure 4 corresponding to v > ^ pea k- We plot the square root of the 

autocorrelation function for the temporal fluctuations of the afterglow flux C( r obs)> as normalized 

1/2 
by the unknown amplitude of the fractional density fluctuations in the ambient medium, £ . We 

normalize tq by Rb, which is different for the three different cases of 7&. 

The value of [C(^obs)] 1 at r b s = provides the typical amplitude of the fluctuations in the 
observed flux. The charateristic period of the fluctuations is 7i/2 j0 bs> the time over which C(^obs) 
drops to half its maximum value. We list the derived values of these parameters in Tables 2 and 
3. The tables implies that: (i) for a given ro/-R&, [C(0)/£o] decreases with decreasing 7^; and 



(ii) for a given 75, [C(0)/£o] decreases with decreasing r^/Rb- Poisson statistics implies that 
density fluctuations with an amplitude <5o <C 1 on a correlation length / will induce an average 
fluctuation amplitude ~ 5q/vN in a region of size L>I, where TV ~ (L/l) s is the number of 
independent regions of positive or negative density fluctuations in the sampled volume. This 
explains the qualitative trend of [C(O)] 1 ' 2 to decrease as the value of r$ is lowered. However, 
K(0)] is n °t proportional to r in our problem, because the different uncorrelated regions 
within the integration volume have different weights in their contribution to the total flux. 



4. CONSTRAINTS FROM GRB990510 

Stanek et al. (1999) have monitored the optical afterglow of GRB990510 on a sub-hour basis 
and obtained an rms scatter of 0.02 mag for the BVRI observations during its first day, consistent 
with the scatter being entirely due to photometric noise (see also Hjorth et al. 1999). These 
observations imply an upper limit of 2% on the rms amplitude of optical flux fluctuations on time 
scales from 0.1 to 2 hours during the first day of the afterglow. After 1.6 days, the afterglow 
decline rate steepened, possibly due to the lateral expansion of a jet (Stanek et al. 1999; Harrison 
et al. 1999). Lateral expansion is expected to be important when the Lorentz factor decelerates 
to a value of order the inverse of the jet opening angle, but could be neglected at earlier times 
(Rhoads 1997, 1999). In the following, we focus on the early stage of this afterglow (< 1.6 days), 
during which the observed afterglow radiation originates from a region much smaller than the jet 
opening angle. At this stage, the observed region behaves as if it is part of a spherically-symmetric 
fireball, and hence should be adequately decribed by our spherical expansion model. Our model 
also assumes that the unperturbed ambient gas has a uniform density. Chevalier & Li (1999) 
studied wind interaction models for GRB afterglows, and concluded that GRB990510 can be 
better explained by a constant density medium than a wind density profile. 

The equivalent 7-ray energy release of GRB990510 for isotropic emission is E = 1.2 x 10 53 erg, 
if the source redshift is z = 1.62 (Wijers et al. 1999). We assume that the hydrodynamic energy 
release is comparable to this value, and also adopt n = 1 cm" 3 . The afterglow emission peaks in 
the optical (f pea k = 5.1 x 10 14 Hz) at T b s = 9.6hr, when 75 = 22 and Rb = 3.9 x 10 17 cm. Since 
most of the observational data from Stanek et al. (1999) is at frequencies v > Vpcak, we compute 
the flux fluctuations in this regime. Our calculations indicate that flux fluctuations on time scales 
( r i/2,obs) between 0.3 and 2 hours correspond to density fluctuations on length scales, ro, between 

3.5 x 10 14 and 2.6 x 10 15 cm, i.e. in the range of ~20-200 AU. Based on the observed upper limit 

1/2 
on the amplitude of flux fluctuations in GRB990510, we calculate the upper limit on £ as a 

function of ro in this range. The resulting constraints on the ISM inhomogenties in the vicinity of 

the progenitor of GRB990510 are illustrated by the solid line in Figure 5. The horizontally-shaded 

region above this curve is forbidden while the region below the curve is allowed, based on the 

afterglow data. We find that the amplitude of density fluctuations has to be lower than ~ 10% 

on the length scale of ~ 200 AU and lower than unity on the ~ 20 AU scale. 
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It is instructive to compare our derived constraints with observational data for density 

fluctuations on similar length scales in the ISM of the Milky- Way galaxy. Structure on small 

scales was first inferred by Dieter, Welch, & Romney (1976), using VLBI observations at 21 cm 

against the extragalctic source 3C 147. Diamond et al. (1989) and Davis, Diamond, & Goss (1996) 

obtained similar results for more sources. Frail et al. (1994) detected temporal variations in the 

21 cm absorption towards six high- velocity pulsars, and inferred changes in the HI column-density 

of ~ 13% on projected scales of 5-100 AU. This sets a lower limit on the rms density contrast 

(5uh I nil) of ~ 13% in spheres of the above scale; the actual density contrast could be much larger 

due to partial cancellations between overdense and underdense regions along the line-of-sight 

to the pulsars. Lauroesch & Meyer (1999) studied the small scale ISM structure in atomic gas 

by observing the interstellar K I absorption line towards multiple star systems, and inferred a 

hydrogen density contrast {8uh jnu) ~ 1 — 2 on the length scales of 10 2 -10 3 AU (Meyer 1999, 

private communication). Small-scale density inhomogenieties were also inferred in molecular 

clouds (Marscher, Moore, & Bania, 1993; Moore & Marscher, 1995). All these studies indicate 

that structure is ubiquitous on scales of 10-10 3 AU in the ISM (for a physical interpretation of the 

above results, see Heiles 1997). The observed structure might be caused in part by fluctuations in 

the ionization fraction or chemistry. However, for the purpose of putting our results in the context 

of these local ISM observations, we will assume that they relate to actual inhomogneities in the 

gas density. In order to compare our results with the above data we need a relation between the 

autocorrelation function and the rms amplitude of density contrast in a region of a given size. For 

a spherical region of radius R, this relation is to a good approximation given by (Padmanabhan 

1993) 

'5n\ 2 \ 1 f R ^ o 9 dk ,,„. 

k 3 4~r , (43) 



\ \ n J r 2tt 2 V Jo 

where V is the normalization volume, a\ is the power spectrum of the density fluctuation and is 
equal to the Fourier transfer of the autocorrelation function: 



a 2 k = V f(r>-* f (?T ? . (44) 

For the autocorrelation function in equation (B2), 

/•OO ("K f'I'K 

a 2 k = V dr dd d^ e- r/ro e- ikrcose r 2 sm6 (45) 

Jo Jo Jo 

= 4ir(,oV I ye~ r/r ° sin krdr. (46) 

Jo k 

Subtituting this result into equation (|43|), we get 



'fcYV^rVr/™*'^ 2 



Co / r e r/r °dr / k sin krdk = —£,o 
n J rI tt Jo Jo vr 



r \ R/r Q 

arctan ' ' 



RJ l + (R/r ) 2 



(47) 



Based on observational constraints for {(6n/n) R ) on a given length scale R, we may now calculate 
from equation (|47|) the related value of £ for any assumed value of tq. 
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The short-dashed curve in Figure 5 shows the constraint imposed by setting {(Sn/n)^) = 1 
for R = 10 2 AU, while the long-dashed curve corresponds to ((dn/n)^} = 1 for R = 10 3 AU. The 
vertically-shaded region between these two curves describes intermediate length scales, refering to 
the range of constraints imposed by the observations of Lauroesch & Meyer (1999). This entire 
region is well above the solid curve calculated from GRB990510, implying that we should not have 
been surprised if flux fluctuations at a level of a few percent were detected in the optical afterglow 
of GRB990510. The dot-dashed curve in this plot corresponds to setting ((^n/n) 2 ^) 1 ' 2 = 0.13 for 
R = 10 2 AU, which refers to the lower limit on the HI density fluctuations inferred by Frail et al. 
(1994). The allowed range of density fluctuations are above this curve. This constraint is close but 
still above the solid curve, i.e. in the region excluded by the variability data on GRB990510. In 
summary, we find that if a GRB occurs in the ISM where the density perturbations have similar 
properties to those inferred by Lauroesch & Meyer (1999) or Frail et al. (1994), then the resulting 
fluctuations in the afterglow flux should be detectable. 



5. CONCLUSIONS 

We have found that linear density fluctuations with 5n/n ^ 1 on the length scale of ~ 1-10 3 
AU could induce afterglow flux fluctuations with a fractional amplitude of up to ~ 40% over 
timescales of tens of seconds in the X-rays, up to ~ 30% over tens of minutes in the optical, and 
up to ~ 9% over days in the radio (see Table 3). These flux fluctuations average over the full 
range of density inhomoegeneities within the emission region. For example, during the optical 
afterglow, the emission region occupies a rather large volume of ~ (10 4 AU) x (10 3 AU) x (10 3 AU) 
(assuming E52 = 1, n\ = 1, z = 1) and so inhomogeneities on scales <C 10 3 AU would surely be 
ensemble averaged as long as their volume filling fraction is not too small. 

At both extremes of high (X-ray) and low (radio) frequencies, the calculated variability 
might be contaminated by other effects. During the early period of the X-ray afterglow, the 
externally-induced fluctuations we considered might be blended with variability associated with 
internal shells within the fireball which are still catching-up with the decelerating blast wave, 
and which were ignored in our analysis. At the opposite extreme of low radio frequencies, the 
flux might scintillate due to inhomogeneities in the local (Galactic) interstellar medium along the 
line-of-sight (Goodman 1998; Frail et al. 1997, 1999; Waxman et al. 1998). In principle, the radio 
flux variability predicted by our model can be distinguished from variability due to scintillations, 
based on its different dependence on photon frequency, especially at high frequencies. However, 
it appears that the best spectral regime to observe the afterglow flux variability predicted in this 
paper is in between the X-ray and radio frequency windows, e.g. in the optical-infrared band. 

Our calculation assumed spherical symmetry and should apply to the early expansion stages 
of a jet, as long as its Lorentz factor is larger than the inverse of its opening angle (Rhoads 
1997, 1999, Sari et al. 1999). Our treatment could be extended in the future to describe the 
lateral expansion of a jet at later times, the possible existence of a power-law density profile as 
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for a precurser wind from the GRB progenitor (Chevalier & Li 1999), and the possible effects of 
nonlinear clumps of density (Dernier & Mitman 1999). 

The application of our simple model to the early optical afterglow of GRB990510 (Stanek 
et al. 1999) provides already interesting upper limits on the density fluctuations on scales of 
~ 20-200 AU around the source. These limits are lower than the observed fluctuation amplitude 
on similar scales in the local interstellar medium (Fig. 5). If these local measurements apply 
to interstellar turbulence in high redshift galaxies, then optical monitoring of future afterglows 
should reveal flux fluctuations at the level of a few percent or higher on the timescale of less than 
an hour. 

We thank Bruce Draine, Dale Frail and David Meyer for useful discussions about the detection 
of density fluctuations in the ISM. This work was supported in part by NASA grants NAG 5-7039 
and NAG 5-7768. 
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lb 


T ohs /[El^n- 1/3 ((l + z)/2)] 


/fc/(j$V /s ) 


*P«*/[nl /2 ((l + *)/2)- 1 ] 


100 


202.6 sec 


6.1 x 10 16 cm 


2.9 x 10 17 Hz 


22 


3.2 hr 


1.7 x 10 17 cm 


6.7 x 10 14 Hz 


3 


27.0 day 


6.3 x 10 17 cm 


2.3 x 10 11 Hz 



Table 1: Values of T b s , Rb and f pe ak for the three choices of jb- ^52, n \ and z are left as free 
parameters. 



case 


7fe 


ro/Rb 


r /(i4 /: V 1/3 ) 


[C(0)/£o] 1/2 


r l/2,obs/^obs 


ri/ 2 ,obs/[^T% 1/3 ((l + ^)/2)] 


1 


100 


0.001 


6.1 x 10 13 cm 


0.05 


0.026 


5.3 sec 


2 


100 


0.01 


6.1 x 10 14 cm 


0.28 


0.30 


61.6 sec 


3 


22 


0.001 


1.7 x 10 14 cm 


0.017 


0.027 


0.086 hr 


4 


22 


0.01 


1.7 x 10 15 cm 


0.22 


0.30 


0.94 hr 


5 


3 


0.001 


6.3 x 10 14 cm 


0.004 


0.026 


0.70 day 


6 


3 


0.01 


6.3 x 10 15 cm 


0.07 


0.26 


7.0 day 



Table 2: Values of [C(0)/£o] 5 r i/2,obs/^obs and Ti/2 i0 bs f° r srx cases in the v < is pea k region. The 
values of E52, n\ and z are left as free parameters. 



case 


lb 


ro/Rb 


[C(0)/£o] 1/2 


r l/2,obs/^obs 


T 1/2 , ohs /[E^n- l/ ' d ((l + z)/2)} 


1 


100 


0.001 


0.06 


0.033 


6.7 sec 


2 


100 


0.01 


0.37 


0.33 


67 sec 


3 


22 


0.001 


0.026 


0.034 


0.11 hr 


4 


22 


0.01 


0.27 


0.31 


1.0 hr 


5 


3 


0.001 


0.009 


0.031 


0.84 day 


6 


3 


0.01 


0.09 


0.24 


6.5 day 



Table 3: Values of [C(0)/£o] ■> r i/2,obs/^obs and Ti/2 i0 bs f° r s i x cases in the u > ^ pea k region. The 
values of -E52, n\ and z are left as free parameters. 
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Fig. 1. — Coordinate system for the integration of the emitted afterglow flux. 




Fig. 2. — Notations and geometry of the integration region. 
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Fig. 3. — The temporal correlation function for variations of the aftreglow flux [C(7"obs)/£o] as a 
function of r b s for v < f pea k> where r b s scales as [E^ % ((1 + z)/2)]. The Six plots correspond 
to six cases for different values of parameters, as shown next to each plot. 
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Fig. 4. — Same as in Figure 3, but for u > v pea k- 



18 



100 



10 



0.1 



0.01 







1 i 


4. 








N 


i 








\ 


I 










I 








\ 


I 














I 












\ 


I 














\ 




















I 


















1 






















M 








































N, 




















N 




















\ 




































""" _ 
































V, 


























^ 












*■» 










x k ^-^- 



\ 



X 



i\ 



Lfiurbesch 
Meyer 1 




Frail et al 
1994 . 



GRB990510 



ALLOWED 



10 1 



10 15 
(cm) 



10" 



Fig. 5. — Constraints on the £o -r o parameter space, for a density autocorrelation function of the 

1/2 

form £(r) = £o ex P( — r / r o)- The solid curve describes the upper limit on £ for r$ between 
4 x 10 14 to 3 x 10 15 cm based on the upper limit on the amplitude of optical flux fluctuations in the 
GRB990510 afterglow (Stanek et al. 1999). The upper horizontaly-shaded region is not allowed 
by the GRB990510 data. The short-dashed curve and the long-dashed curve reflect constraints on 
the inhomogeneities in the local ISM, based on Lauroesch & Meyer (1999), with the short-dashed 
curve corresponding to {(5n/n) 2 R ) = 1 in spheres of radius R = 10 2 AU, and the long-dashed curve 
corresponding to ((5n/n) 2 R ) = 1 in spheres of radius R = 10 3 AU. The vertically-shaded region 
between these two curves describes intermediate length scales. The dot-dashed curve corresponds 



to {{dn/n) 2 ^) 1 ' 2 = 0.13 in spheres of radius R 
inferred from Frail et al. (1994). 



10 AU and describes the lower limit on £ 
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